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Abstract 

This article will review recent results on dimensional reduction for 
branched polymers, and discuss implications for critical phenomena. Parisi 
and Sourlas argued in IPS81I that branched polymers fall into the universal- 
ity class of the Yang-Lee edge in two fewer dimensions. Brydges and I have 
proven in |BIfll| that the generating function for self-avoiding branched poly- 
mers in D -I- 2 continuum dimensions is proportional to the pressure of the 
hard-core continuum gas at negative activity in D dimensions (which is in 
the Yang-Lee or itp'^ class). I will describe how this equivalence arises from 
an underlying supersymmetry of the branched polymer model. 

I will also use dimensional reduction to analyze the crossover of two- 
dimensional branched polymers to their mean-field limit, and to show that 
the scaling is given by an Airy function (the same as in ICarOlI 'l. 

Let us begin with a definition of branched polymers, first on the lattice, 
then in continuous space. On the lattice, branched polymers should look like the 
following picture. 



Define cn to be the number of iV-vertex branched polymers mod translations. 
Let r be a tree graph on {1, . . . , N} and yi, . . . ,yN be an embedding of the graph 
in the lattice. In other words, yi denotes the position of the i*"^ vertex of the tree 
graph. If i and j are connected by a line in the tree graph, then yi and yj must be 
adjacent in the lattice. In addition, no two vertices may be mapped to the same 
point; this would lead to a loop. The generating function Zbp{z) = J2n ^^^n can 
be written as 

N=l T V2.---,yN ij£T ij^T 

where 

2[/^ =(5|j,,_j,^l,i (2) 

and 

Utj = 1 - 5y^,y, (3) 

enforce the adjacency and loop-free conditions, respectively. For example, C3 = 
^ • 3 • 4 • 3 = 6 in as there are 3 tree graphs on {1,2,3}, and for each of them 
there are 4 choices for y2 and 3 for y^. See |BI01I llYoSG) for details. 

Background 

The concept of dimensional reduction was introduced by Parisi and Sourlas in 1979 
in order to understand the behavior of the Ising model in a random magnetic field 
|PS79) . Although the proof of long-range order in three dimensions |Imb84l Iimb85| 
showed that dimensional reduction did not work in that problem, dimensional 
reduction appeared to work for branched polymers |PS81| . Our proof answers any 
lingering questions about dimensional reduction for branched polymers and gives 
a clear picture of how it works. The question of what goes wrong for random fields 
continues to be debated in the literature |BD98llP^02llFel02| . 

Dimensional reduction for branched polymers was employed recently by 
Cardy to determine the crossover scaling function for area-weighted self-avoiding 
loops in two dimensions |Car01) . Below, I will use our dimensional reduction iden- 
tity to derive a similar result for the crossover from branched polymers to their 
no ninter acting or mean-field limit. 

Branched Polymers in M^+^ 

Lattice models of branched polymers have been written down that are supersym- 
metric after dropping irrelevant terms |Sha83irSha85| . Here we work in continuous 
space so that the supersymmetry responsible for dimensional reduction can be 
preserved. Let us adjust the definitions given above for lattice branched polymers 
to the continuous case. In R^+^, a branched polymer consists of 

1. A tree graph T on 1, . . . , TV and 
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2. An embedding y = {yi}i=i^....N into R^+^ such that 

• if ij G T then \yij \ — 1 and 

• if ij ^ T then \yij \ > 1. 

The generating function for branched polymers (mod translations) is 



^Bp(z) = E In E / • • • dyiv n n f^^^' (4) 

where 

U^,=U{\y,-y,\') (5) 

and 

2C/4.=2(7'(|y.-2/,f) = 5(|y.-2/,|-l). (6) 

As before, Ul^ enforces adjacency, and Uij is a self-avoidance condition. One can 
think of these branched polymers as collections of hard spheres with kissing con- 
ditions determined by the tree graph T. 

The Hard-Core Gas in 

We shall relate branched polymers in R^"'"^ to the hard-core gas in MP . The hard- 
core gas has hard spheres centered at Xi € MP for i = 1, . . . , A^. Let 

and let 

u(kj) = e{Uj - 1) (8) 

enforce the hard-core constraint. Then 

^HcW=E^/ dx^---^^f^ n V{\x.,-x,f). (9) 

Main Results 

Theorem 1. For all z such that the right-hand side converges absolutely, 

lim -LiogZuciz) = -27rZBp(-^) . (10) 
A/Ro |A| V 27r/ 
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Examples 

Exact computation of the left-hand side is possible in D = 0,1. For D = 0, 
Zhc{z) = 1 + 2 and so in dimension D + 2 = 2, 

^BpW = -^log(l-27rz). (11) 

The pressure of the one-dimensional gas of hard rods is also computable, see for 
example |BI03| . Hence in dimension D + 2 = 3, 

Zbp(z) = -T{2nz) = Y: ^(2.z)^. (12) 

N=l 

Here T{z) ~ — LambertW(— z) is the tree function jCGHJK] , These equations 
allow us to write a simple expression for cjv in two and three dimensions. 

Green's Functions 

We also have a dimensional reduction formula for correlations. First, some defini- 
tions. Let 

N N 

Pix) ^YS{x - x^), p{y)^YS{y~y^), (13) 

i=l i=l 

where i, Xi G and y, yt G Then the density-density correlation functions 

of the two systems can be written as 

GBp{yi,y2;z) = ^ ^ / dyi • • • dyjv p(yi)p(y2) H t^^y] IT ^y'' 

N=l ■ T •' ijeT ij(^T 

Ghc{xi,X2;z) = lim (p(ii)/)(i2))Hc A ■ (14) 

Here (•)hc,a is the expectation in the measure for which Z-aciz) is the normalizing 
constant. 

Theorem 2. 

Ghc(0, X- z) = -2tt j (fwGBP (O, y; -|^) , (15) 

where y — (x, w) <E K'°+^. 

Analogous results can be obtained for the n-point functions. 
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Consequences for Critical Exponents 

One can define a critical exponent ajjc by looking at the singularity of the pressure 
of the hard-core gas: 

Likewise, a susceptibility exponent 7hc can be defined from the behavior of Zbp : 



^BP 



-^)-{z-z,)^-''--. (17) 



Keep in mind that Zc is negative and — is positive — Theorems ^ and |21 relate 
the critical behavior of branched polymers to the critical behavior of the hard-core 
gas at negative activity. Theorem ^ equates the two singularities, so 

OHC = 7BP- (18) 

Examining the 2-point functions of each model near their critical points, we 
can define exponents ?7hc7 ^bp from the asymptotic forms 

Ghc(0,x;2) ~ x-(^-2+„Hc)j^^jj^(^/^jj^)^ (19) 

Gbp(0,x;z) ~ x-('^-2+"Bp)A3p(a;/^3p) (20) 

as \x\ ~* oo with x := x/^ fixed. Here d = D + 2, ^hc, Cbp arc the correlation 
lengths of the respective systems, and Kuc, Kbp are scaling functions. Theorem 
121 implies that ^hc = ^bp, so in particular the exponents vnc = i^bp, where 

eHc(^) - {z^zc)-^-, (21) 
Cbp(z) ~ (z-ze)-''-- (22) 

as z \ Zc. Theorem 13 also implies that 

mc = VBP (23) 

for D > 1. Furthermore, if one defines 9 from 



JV=1 

with 



Zbp(z) = ^^^"^ (24) 
CN ~ z-^N-''. (25) 



Then 

= 3- 7BP. (26) 
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Further results can be obtained for the scahng functions. We can express 
J^BP in terms of -fCnc ^ind its radial derivative: 

KBp{i) = ^ [ii^HcC^) -{D-2 + ?/Hc)i^Hc(i)] • (27) 
In £) = 1 we can compute |BI03) 

Kkc{x) = -^ e-^ (28) 

so that in + 2 = 3 dimensions, 

KBp{i) = e-\ (29) 

which agrees with the prediction of Miller |Mil91| . Let us tabulate the known or 
for the critical exponents of the hard-sphere gas in various dimensions. 



D 


a 




77HC 


a = 1 — a 





2 






-1 


1 


3 
2 


1 

2 


-1 


1 

2 


2 


7 

6 


5 

12 


4 

5 


1 

6 


MFT D>Q 


1 
2 


1 
4 





1 
2 



In £) = 0,1, the values are rigorous. In D = 2, the value a = ^ is obtained 
from Baxter's solution to the hard hexagon model |Bax82l IDha83l lBL87| . The 
last column gives values for cr, the Yang-Lee edge exponent. Assuming that the 
hard-sphere gas at negative activity falls into the Yang-Lee edge universality class, 
we have a = 1 — a. Actually, from hyperscaling (Di^hc = 2 — anc) and Fisher's 
relation know that there is only one independent exponent in 

each dimension. The value rync = ~ f independently confirmed from conformal 
field theory |Car85j . 

The branched polymer critical exponents in D + 2 dimensions can be equated 
to the corresponding values for the hard sphere gas in D dimensions; see the first 
three columns of the following table. 



D + 2 


7 = a 




?7BP = '7HC 


61 = 3-7 


2 


2 






1 


3 


3 
2 


1 
2 


-1 


3 
2 


4 


7 

6 


5 
12 


4 
5 


11 

6 


MFT D > 8 


1 
2 


1 
4 





5 
2 
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The values in two and three dimensions are now rigorous. Note that in high di- 
mensions {d > 8) it has been proven that 7bp = 5 , fsp = j , Vbp = (at least for 
spread-out lattice models) |HS90I I11S92I IHvS03j . While our results do not apply 
to lattice models, they give a strong indication that the corresponding hard-core 
exponents have the same (mean- field) values for D > 6. 

Forest-Root Formula 

The underlying mechanism behind all these results is an interpolation formula 
(the Forest-Root formula). Let /(t) depend on t = {tij), (ti) for I < i < j < N. 
Assume that / at 00. Let = \wi — Wj\'^, ti = \wi\'^ with Wi G C Then 




The sum is over forests F and roots R (collections of bonds ij, and vertices i, 
respectively) such that each tree of F has exactly one root R. See Fig 1. 




tree 



tree 



tree 




tree 



forest 

Figure 1: Example of a forest 

For a simple example, consider the case N — 1,F = 0,i? = {1}. Then the 
formula reduces to the fundamental theorem of calculus in one variable: 



/(0)= / f'{t)^ = - f'{t)dt. 



(31) 



Supersymmetry 



To understand the Forest-Root formula in general, we need to exploit supersym- 
metry. Replace in /(t) each variable ti with 



Ti = WiWi 



dwi A dwi 



(32) 
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and tij with 



(Recall that Wi. 



dwij A dwij 



(33) 



Wi 



-Wj). Then /(r) is defined by its Taylor series. 1 claim that 



a "localization" formula holds: 



/fe) = /(o). 



(34) 



This formula becomes the Forest-Root formula when expanded out. The absence 
of loops comes from the fact that {dw A dw)^ = if G has a loop. We prove 
the formula by deforming the problem to the independent case, using ideas from 

jWltSU- 

An alternate argument can be given, exploiting the linearity of (|34|) to reduce 
to the case where / is an exponential |BW88| . Let 



(35) 



/(a) = / dtf{t)e'' 



be the Laplace transform of /. Here a = (oi), (a^ ) is dual to t, that is a • t = 
TiijQijtij + l^iQiti. Let us assume that / has exponential decay in each ti variable, 
so that we can take the a^'s to have positive real part. Cases where / decreases 
more slowly at infinity can be treated by taking limits. We can write / in terms 
of its Laplace transform. 



fit) 



zm zni 



(36) 



with the integration contours running parallel to the imaginary axis. We see that 
H34|l holds if it is true for functions of the form e^"* *. 
Let us write 



a • t = (w. Aw) 



dw , dw 
,A- 



'2m 



where 



Aij — Qij , 



(37) 

(38) 
(39) 



The matrix A is symmetric and has positive real part. Then by the properties of 
the differential forms dwi, dwi, 



^— (w^Aw) 



- {w,Aw) 



exp 



dw . dw 
--,A- 



N 



detAjl 



i=l 



—dwi A dW; 
27ri 



(40) 
(41) 
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With Wi = Xi + iyi, we have that dwt A dwi = ~2idxi A dyi, so this is 

^-^det A J e-^-^^^^^d^xd^y = 1, (42) 
and we obtain for /(t) = e^"* *. 
Decoupling in Two Extra Dimensions 

Use the Forest-Root formula to decouple the spheres in the hard-core gas by mov- 
ing them apart in two extra dimensions. At fixed N we have an integral over 
X € C M^^ of 

/(0)= n Ui\x,,\'), (43) 

l<i<j<Ar 

where ti = |u'ip = 0. Extend this to it) 7^ by writing 

/(t) = n Ui\x,,\^ + Uj) X (large t cutoff). (44) 

l<i<j<A' 

Apply the Forest-Root formula, 




Then each , when applied to Uij , becomes ^ surface measure for the combined 
integrals over = (xij , Wij ) . This means that spheres are stuck together accord- 
ing to the forest F. Furthermore, the trees of the forest decouple in the limit as 
the large t cutoff is removed. This is because they are spread apart in the w di- 
rections, so that they rarely interact. As a result, Zhc can be written in the form 
of an exponential of tree graphs. Thus log Zhc is evaluated as a sum/integral over 
branched polymers, and Theorem ^ is proven. See |BI01| for details. 



Relation with the Yang-Lee Edge 

Repulsive gases near criticality at negative activity are described by an iip^ field 
theory. More precisely, they should fall into the same universality class — see jCar82l 
IDha83l ILF95I IFFOO] . To see this, replace the hard-core potential with a repulsive, 
smooth, rapidly decaying two-body potential 



u{xij) = t'(|xyp) with u{k) > 0. 



(46) 



By the sine-Gordon transformation, the repulsive gas 



Zy{z) 



y — 

JV=0 



dx\ ■ ■ ■ dxN exp 



-/? ^ U{X^J] 
l<i<j<N 



(47) 
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can be written as a — ze*''' field theory. With z = ze^'"''^^/'^ , we have 

1 / s-l \ 

Zv{z) = y "sxp 

If z is negative enough, this action looks critical. In fact at the critical value 
of z, the lowest order term in the action is iLp^ . A similar situation occurs for Ising 
or ip^ model in a sufficiently large imaginary field |Fis78| . 



e 2 



(48) 



Dimensional Reduction: Soft Branched Polymers 

The arguments leading to Theorems ^ and |3 work for any U{t) for which U — \ 
and U' have sufficient decay. Take 

U{t) = exp(-/??;(t)) (49) 

to get the identity 

where Z^y',v is given by the same formula as before but with this U . In this 
model of "soft" branched polymers, the repulsion between monomers is given by 
a factor exp(— /3t;(|xij p)), and the attraction linking neighboring monomers in 
the tree graph is given by a factor —2f3v'{\xij\^)). Together with the sine-Gordon 
representation for Zbp^v, we obtain a direct representation for branched polymers 
in £'-1-2 dimensions as a —ze'^'^ field theory in D dimensions. As this is presumably 
in the Yang-Lee class, we have confirmed the prediction of Parisi and Sourlas 

lEHHi]. 



Crossover to Mean-Field Branched Polymers 



We wish to understand what happens as the repulsion between monomers is turned 
off. In the limit, we obtain a mean- field or noninteracting model (on the lattice 
such models were analyzed in jFro86llBCvSj ^. Let v = v[Q) be small and z be large 



so that a redefined z 



zve 



v/2 



is fixed. In these variables (|50|l becomes 



log 



exp 



V 



dtp 



~ — 27rZRP.„ I — 

ZTT 



(51) 



I have reversed the sign of z, so the branched polymer singularity is now at positive 
z. 

It should be clear that the w ^ limit simplifies both sides of this identity. 
On the 0-dimensional Yukawa gas side, w — *■ means the integral reduces to finding 
the critical point of On the two-dimensional branched polymer side, 

V ^ Q eliminates the self-avoidance. Every tree has one fewer edges than vertices. 
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Assuming v' remains proportional to w, as w — *■ and zv — > z, this means that 
j^Zbp^v goes to a finite limit, which is a generating function of non-self-avoiding 
(or mean-field) branched polymers containing the origin. 

We can observe the crossover by considering function of z and v 

near v = and z = z^- With S{lp) = ^tp^ -\- ze'"^, we have 

S'{(p) — ip + ize^^ — when iip ~ x with xe~^ ~ z (52) 




The solution is a; = x{z) — T{z) as in the graph. We also solve 

S"{ip) = 1 - ze^'^ = when z = e-'^ = e"^^ = : 



(53) 



so = 1, 5c = e-i has S" = S" = 0. 

Let us expand S{lp) about ip = —i with z close to Zc = e^^. With := (/? -I- i 
we obtain 



Let t 



'-. This becomes 



k + ez e'^ 



(p^ -i(p^\- e(zc - z)e'^ + e 



1,-2 
2 



(54) 



(e''^ - 1 - ip^\pP-)-t{e'^-\)-t^\ = -'t^-itp^O{ip^)+tO{p?)-t^\, (55) 
and so from H51(l 



-27rZBP,„ ( ^ ) = log / '^^P 



o3 



i ( ^ + + 0{p^) + tO(^^) + 1 - ^ 
w V 6 2 



/27rw 
(56) 

Let us apply ^ to get the sum of branched polymers containing the origin. With 
G^7..v{y) = J;2'Bp,t,(?/), we have 



c?z 



log / exp 



1 / ip' 



:;3 



f V 6 



-i-ii^ + 0(^*)+iO((^^) + t 



(57) 



The expression in the exponent is useful for p small. Consider the crossover 
regime, \t\ < 0(w^/'^). I wish to show that the contribution to the integral from 
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\(p\ > w^/^""^ is smaller than any power of v. Make a kink in the contour of inte- 
gration for (f, so that it lies at angles |- and This turns into — From 
H55|l . one can check that on this contour, 



ReS{ip) > const min{|(^p, - tOifp) - t + 5 



(58) 



so that 1^ is effectively of order v^^"^ . This means that ^Oiip^) and ^0(<p^) are 
0[v^l^^ and that these terms do not contribute to Zbp.v to leading order in v^^^. 
Likewise when computing ~ ve^^^j^ = —eve"^^-^ in H57|l . the term 0{(f'^) = 
0{v'^/'^) is subleading order. The term — ^ ^ = e^+^Z^ is nonsingular. Hence the 
singular behavior of Gbp,^ is given to leading order by 



^(sing) / Z 



dz 



log / exp 



~6~ 



itip 



dip. 



(59) 



We have arrived at a point very similar to the one Cardy reached in his 
analysis of the crossover from area-weighted self-avoiding loops to self-avoiding 
loops |Car01| . He related that problem to branched polymers and, as wc have 
done here, used dimensional reduction to compute the crossover scaling function 
as an Airy integral. If we rcscale ip — > (pv~^^^ in H59(l . we obtain as in |(]ar01| that 



G 



(sing) 
BP.v 



v^F{tv 3), 



where 



ds 



(60) 



(61) 



Here 60 and 61 are constants, and Ai{t) = ^ e^^'f^^'f^ I'^dip is the Airy function. 

The expression (|60|l interpolates between two different critical behaviors for 
branched polymers. We have already seen that 7bp = 1 in two dimensions; this 
corresponds to the pole in (|60|) that occurs when v is fixed and nonzero, and t 
decreases through to the first zero of P^\(\)\tv~^\ This is the limit of convergence 
of .Z'bp; past this point H51|l no longer makes sense. 

If, instead, we send v to zero and then let t \ 0, this means Cp is replaced by 
its critical value (pc where S'{(pc — «) = 0, c.f. H52|l . Then, 



lim Gbp t, f— 1 = Um — ue'"/^4- log / exp 



V 



dp 



(62) 



It is clear from the graph above that iipc — t^/^. Thus the leading singularity 
of ()62fl is ' — t^^^, and ^ = \- T his agrees with the value determined for the 
mean field models in |Fro86llBU^ . This also agrees with the conjectured value 
of the entropic exponent of self-avoiding walks. The square-root behavior can be 
seen in (|Hn|) also, as F{s) ~ s^^^ as s 00. 
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The identity of the crossover scahng functions and critical exponents for self- 
avoiding loops and branched polymers suggests a close connection between the 
two systems. It is natural to think of loops arising as the frontier of branched 
polymers. In this picture, the sclf-avoidance interaction of the branched polymer 
gives rise to an effective area weighting for the loop, and the mean-field limit of 
branched polymers corresponds to loops without area weighting. This fits in with 
recent results |LSW01| on the Hausdorff dimension of the frontier of Brownian 
motion (suggesting a link with self-avoiding walks). 
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